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ABSTRACT 

Sightlines through high- and intermediate-velocity clouds allow measurements of 
ionic gas phase abundances, A, at very low values of H I column density, N(H I). 
Present observations cover over 4 orders of magnitude in N(H I). Remarkably, for 
several ions we find that the A vs N(H I) relation is the same at high and low column 
density and that the abundances have a relatively low dispersion (factors of 2-3) at 
any particular N(H I). Halo gas tends to have slightly higher values of A than disk gas 
at the same N(H I), suggesting that part of the dispersion may be attributed to the 
environment. We note that the dispersion is largest for Na I; using Na I as a predictor 
of N(H I) can lead to large errors. 

Important implications of the low dispersions regarding the physical nature of 
the ISM are: (a) because of clumping, over sufficiently long pathlengths N(H I) is a 
reasonable measure of the local density of most of the H atoms along the sight line; 
(b) the destruction of grains does not mainly take place in catastrophic events such as 
strong shocks, but is a continuous function of the mean density; (c) the cycling of the 
ions becoming attached to grains and being detached must be rapid, and the two rates 
must be roughly equal under a wide variety of conditions; (d) in gas that has a low 
average density the attachment should occur within denser concentrations. 

Subject headings: ISM: clouds, ISM, Abundances Galaxy: halo, Galaxy: disk 

1. Introduction 

The study of elemental abundances in the interstellar medium has a long history, both in 
the optical and with space-based UV spectrographs, such as Copernicus and the Hubble Space 
Telescope (HST). These studies show that most elements have apparent abundances relative to 
hydrogen substantially below those in the Sun. This is interpreted as due to depletion of ions onto 
dust grains. 

Savage & Sembach (1996a) summarized high-quality abundance studies based on HST data. 
They concluded that the depletion of an element depends on the nature of the interstellar medium 
along the sightline. In cold, dense disk gas, S has near-solar abundance, while Si, Mg, Mn, Cr, 
Fe and Ni are progressively more depleted. In warm-disk gas, the gas-phase abundances appear 
higher than in the cool gas, while in halo gas they are even larger. 
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Jenkins, Savage & Spitzer (1986) used Copernicus to measure accurate column densities 
for Mg II, P II, CI I, CI II, Mn II and Fe II, for stars with distances, d, between 0.1 and 2kpc. 
They measured N(H I) from the damping wings of Lya and found a relation between the ionic 
abundance relative to H I, A, and the average density, («h) (= N(H l)/d). An explanation was 
provided by Spitzer (1985), who suggested that the observed depletion in a line of sight depends 
on the fraction of warm and cold (dense) gas encountered, where each component has a specific 
gas phase abundance of each ion. 

Crinklaw, Federman, &: Joseph (1994) observed Ca II and Ti II toward 12 stars with distances 
between 150 and 650 pc. The slope of the relation between log[(nn)] and log[A(ion)] is similar for 
Ca II and Ti II, but it differs from those for Fe II, Mn II, Mg II and P II. They proposed that Ca II 
and Ti II mostly occur in the warm intercloud medium, and that the apparent depletion is very 
sensitive to the inclusion of high-density regions (with high depletion) in the sightline. 

These studies were limited to sightlines with N(H I) in the range 10 20 to 7xl0 21 cm~ 2 . To 
extend this range, we used the high- and intermediate-velocity gas (HVCs and IVCs, see Kuntz &: 
Danly 1996, Wakker & van Woerden 1997, Wakker 2000). In 21-cm emission these stand out in 
velocity from gas in the Milky Way disk, so that weak components (down to 10 18 cm -2 ) can be 
measured. At low velocities, parcels of gas with such low column densities can not be separated 
from other high column density concentrations in the line of sight. 

Wakker (2000) analyzed all the published absorption-line data pertaining to the high- and 
intermediate-velocity gas. Improved H I column densities (and thus improved ion abundances) 
were determined for about 250 sightlines. Ca II was measured for 97 components, Ti II for 6, 
Fe II for 13, Mn II for 5, Mg II for 10, and Na I for 34 (low-velocity Na I was taken from Ferlet 
et al. 1985). Most IVCs have intrinsically near-solar abundance (Wakker 2000). However, HVC 
complex C has yl~0.1 solar (Wakker et al. 1999), and the Magellanic Stream has A~0.25 solar 
(Gibson et al. 2000). Including these HVC complexes could lower some of the apparent abundances 
(in 8 cases for Ca II, 4 for Fe II and 10 for Mg II), but no clear effect is visible. [An anomalously 
low value for Mg II, in the VHVC probed by Mrk205, was excluded as this cloud may have low 
intrinsic abundance, and N(H I) is uncertain.] 

No distances are known for many of the HVCs/IVCs, and they fill only a small fraction 
of the pathlength. Thus, we looked at the relation between A and N(H I), rather than («h)- 
Unexpectedly, we found tight anticorrelations of N(H I) with the abundances of Mg II, Ca II, 
Ti II, Mn II and Fe II. A minimal correlation is expected because both quantities are integrated 
along lines of sight with varying physical conditions. Variations include (a) the ionization of 
H, (b) widely differing histories of grain destruction by interstellar shocks, and (c) considerable 
differences of the interstellar radiation field impinging on the gas. Ca II and Na I are not even the 
dominant stage of ionization in neutral regions, so the radiation field is of critical importance for 
their ionization fraction. In §2 we show the correlations between column density and abundances, 
while in §3 we discuss some implications. 
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2. Correlations between abundances and N(H I) 
2.1. The data 

Figure 1 shows scatter plots of log[N(H I)] vs logL4(ion)], the gas-phase abundance relative to 
H I, for Mg II, Ca II, Ti II, Mn II, Fe II and Na I. The least-squares fit and dispersion (labeled 
"rms") are given in each panel. Symbols indicate the origin of the datapoints, as detailed in the 
figure caption. The Local ISM points (open squares) were not used in the fits (see §2.5 below). 
The high-quality data obtained with the "Goddard High Resolution Spectrograph" ( GHRS) are 
presented separately, in order to show the difference between halo and disk sightlines more clearly 
(§2.5). 

The currently preferred oscillator strengths of Mn 11-1197, 1201 and Mg 11-1239, 1240 are 0.20 
and 0.67 dex higher than those used by Jenkins et al. (1986) (see Savage &: Sembach 1996a). We 
therefore corrected their Mn II and Mg II column densities downward by these amounts. 

As the highest abundances are associated with high-velocity gas, we might be seeing 
a manifestation of the Routly & Spitzer (1952) effect (this is the increase in the ratio 
N(Ca Il)/N(Na I) with LSR velocity). This effect is interpreted as showing that Ca is less depleted 
at higher peculiar velocities. However, it is based on nearby (<100pc), low-velocity (<20kms _1 ) 
gas (e.g. Vallerga et al. 1993). For the HVCs/IVCs the LSR velocity is not a good measure of the 
peculiar velocity relative to their surroundings. Further, the IVCs for which both Na I and Ca II 
have been measured do not show the effect (Wakker 2000). 

We now discuss some implications of the i vs N(H I) relations. 

2.2. Result 1: small dispersions 

For Mg II, Mn II and Fe II the standard deviations of logL4(ion)] from the least-squares fit 
at a given N(H I) are only 0.27 in the log. That is less than a factor 2 either way, even though 
the correlation extends over 4 orders of magnitude in N(H I) and over 2 orders of magnitude in 
^4(Fe II). The standard deviations are slightly larger for Ca II and Ti II (0.4 in the log, or a factor 
2.6 either way). The largest scatter is seen for Na I, where it is 0.52 in the log or a factor 3.5 either 
way. 

The larger scatter for Ca II and especially Na I may be related to the fact that, unlike the 
other ions, these are not the dominant ionization stage in the diffuse ISM, so that ionization effects 
may play a larger role. However, Ti II is the dominant ionization stage, and it shows a scatter 
that is comparable to that of Ca II. 

Ferlet et al. (1985) proposed a slope near zero for the relation between Na I abundance and 
N(H I) for the low- velocity gas. This has been used extensively to estimate N(H I) from N(Na I). 
The lower column density IVC data tend to show slightly higher abundances. However, the spread 
in A(Na I) at any given value of N(H I) is rather large. For a given N(Na I), the derived N(H I) 
will be a factor >3.5 off either way in 33% of the cases. The most deviant points in the Na I 
diagram differ from the mean relation by a factor >20. Thus, although on average N(Na I)/N(H I) 
is independent of N(H I), N(Na I) is a rather poor predictor of N(H I). 
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2.3. Result 2: differing slopes 

The slopes of the correlations differ substantially: -0.78±0.04 for Ca II, -0.69±0.08 for Ti II, 
-0.59±0.04 for Fe II, -0.39±0.04 for Mn II, -0.24±0.04 for Mg II, and -0.16±0.06 for Na I. 
We find slopes of — 0.63±0.12 and — 0.66±0.23 for low- and high- velocity Ca II; for Fe II we find 
— 0.54±0.05 and — 0.70±0.07, respectively. For other elements there are insufficient datapoints to 
make separate fits. The errors are larger because of the reduced range in N(H I), but the resulting 
slopes are the same to within the formal error. This shows that on average the high column 
density, low-velocity disk gas behaves in the same manner as the low column density high-velocity 
halo gas. 

2.4. Result 3: no obvious ionization effects 

The ionization potentials (LP.) of Fe II, Mn II, and Mg II are 16.2, 15.6, and 15.0 eV, 
respectively (compared to 13.6 eV for H I), so these ions can co-exist with H + . The Fe IIl/Fe II 
ratio should not depend on N(H I), but only on the surrounding radiation field (which only 
depends on location). Thus, if the interstellar radiation field has similar values in the sightlines 
that were observed, and if H + were only present in shells around H I cores, then the fraction of 
H + should increase with decreasing N(H I), resulting in an overestimate of the abundance, when 
calculated as N(ion)/N(H I). This would produce an upturn at low N(H I) in the A vs N(H I) 
relation. However, such an effect is not obvious. Ca II (LP. 11.9 eV) should ionize more easily than 
H, leading to a downturn at low N(H I). This is also not seen. 

2.5. Result 4: environmental dependence 

For Mg II, Mn II and Fe II the GHRS disk points (open circles in Fig. 1) scatter around 
those of Jenkins et al. (1986). GHRS halo sightlines (filled circles) tend to lie above the average 
correlation. The small scatter in ^4(Fe II) remarked upon by Savage & Sembach (1996a) is 
consistent with the fact that their halo points span only a factor 10 in N(H I). However, some 
more recent halo points (filled triangles) have higher Fe Il/H I ratios, as do some HVCs. 

The low N(H I) Local ISM (LISM; open squares) shows abundances far below halo points at 
the same N(H I), although they are near the high end for disk gas and the observed depletion is 
consistent with the average Local ISM density of 0.1 cm" 3 . We surmise that both the Sun and 
the LISM probes lie within a cloud, so that the LISM paths do not pass through all of it, unlike 
the case for clouds in the sightlines to more distant stars. Thus, the Local ISM points suggest 
that density is the relevant physical parameter determining depletion. Over sufficiently long 
pathlengths column density can substitute because it is dominated by the densest regions. 

The largest depletions are seen toward ( Oph and £ Per, which were used by Sembach & 
Savage (1996) as the archetypical cold disk sightlines. These have high N(H I) and A(Mg II) and 
A(Fe II) lie below the average relation. Including N(H2) would move these points 0.2 and 0.4 dex 
to the right and down, respectively, shifting them ~0.1 and 0.2 dex further below the average. 
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We thus see the environmental dependence found by Sembach & Savage (1996): disk sightlines 
tend to have low gas-phase abundances, while halo sightlines show high gas-phase abundances. 
We can understand this in large part as a column density dependence, as disk sightlines tend to 
have higher column density than halo sightlines. 

3. Discussion 

The general trend that the gas phase abundances of refractory ions decreases with N(H I) is, 
of course, well known; the ions are removed from the gas phase by sticking to grains. However, 
the low dispersions of the depletions at a particular value of N(Hl), continuing to low values to 
N(Hl), are remarkable. The depletion patterns provide a diagnostic of the density fluctuation 
within the ISM, since the depletions require encounters between the grains and ions. The low 
dispersions of A at a given N(Hl) suggest additional properties of the ISM: 

1. Large values of N(H I) arise from one or a few regions of relatively large local density 
rather than from a superposition of many regions of lower density. In this way, there is unexpected 
physical meaning to the integrated column density: a good correlation with the mean local density 
of the typical H atom along the sightline. If large values of N(H I) were built up from many regions 
of low density along the line of sight, the averaged A would correspond to the lower densities. 
The association of large N(H I) with a large local density is strengthened by the fact that the 
anticorrelation of gas phase abundances with (rin) is better than with N(H I) (Savage & Bohlin 
1979). 

2. The low dispersions in logL4] suggest that grains are not destroyed in strong shocks, as 
some models have suggested. If such destruction occurred, a few points with large gas phase 
abundances would be expected at large N(H I). Rather, the detachment of a refractory element 
from grains seems to be a continuous function of the mean density, and is common enough to 
produce a low dispersion about the mean. Grain destruction may occur by means of many small 
shocks with a broad distribution of velocities. 

3. Both the rates of ions sticking to grains and becoming detached from them (probably in 
atomic form) must be rapid enough to achieve rather good statistical equilibrium. These balancing 
rates lead to a particular value of A for a given density; otherwise, the past history of the H I 
along the line of sight would strongly influence A. 

4. In view of (3) above, we think that the best possibility is that attachment takes place 
within regions of far greater than average density distributed along each line of sight. Encounters 
between grains and ions are very slow if the mean density of the diffuse ISM is assumed (see 
below) . Possibly magnetohydrodynamic waves create non-thermal motions of the ions relative to 
the grains, since the charge per mass of the grains is vastly less that for the ions. As a result, ions 
encounter more dust grains and the net attachment rate increases. 

Under the assumption that relative motions were thermal, Weingartner & Draine (1999) 
calculated average gas phase abundances of a highly depleted species (Ti II) for large N(H I). In 
their model, dense ((nn)~10 3 cm -3 ) molecular clouds contain half the mass of the ISM, 30% is in 
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cool, neutral clouds ((nn}=30 cm~ 3 ) and 20% is in warm neutral clouds ((nH)=0.4cm -3 ). The 
attachment of the Ti onto grains occurs in the molecular and cool clouds. This model shows that 
the observed A(Ti II) for large N(H I) can only be achieved by circulation into very dense regions. 
To achieve rapid attachment within the diffuse ISM, such circulation seems necessary, operating 
at both low and high column densities. An additional condition for rapid attachment of ions is 
that the gas be cold as well as dense, in order to provide efficient focusing by the Coulomb force 
onto negatively charged grains. 

We appreciate comments and clarifications by an anonymous referee, Alex Lazarian, Chris 
Howk, and Bruce Draine. B.P.W. is partially supported by NASA grant NAG5-9179. 
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Figure 1. Correlation plots of the abundances of Na I, Mg II, Ca II, Ti II, Mn II and Fe II vs 
N(H I). The solid lines show the least-squares fits, whose coefficients are given by the label in the 
lower left corner, "rms" is the rms of the residual, also indicated by the two parallel dotted lines, 
"p" is the correlation coefficient. Filled circles are for the HVC and IVC data from Wakker (2000). 
If measurement errors were given for N(ion), error bars are shown; usually no errors are available. 
Crosses show data for low- velocity gas from Ferlet et al. (1985; Na I), Jenkins et al. (1986; Mg II, 
Mn II and Fe II) and Crinklaw et al. (1994; Ca II, Ti II). Errors are usually smaller than the size of 
the dot. The right panels for Mg II, Mn II and Fe II separately show measurements with relatively 
low systematic errors obtained with the GHRShy Savage et al. (1993), Sofia et al. (1994), Cardelli 
et al. (1995), Spitzer & Fitzpatrick (1995), Savage & Sembach (1996b) and Sembach & Savage 
(1996). For these data, halo sightlines are shown by filled circles, disk sightlines by open circles. 
Filled/open triangles show halo/disk sightlines from Fitzpatrick & Spitzer (1997; after combining 
their fit components into the 5 H I components visible in the spectrum) and Howk et al. (1999). 
Open squares show Local ISM sightlines from Linsky et al. (1995), Hebrard et al. (1999) and 
Howk (priv. comm.); these were not used in the fits, see §2.4. For all low- velocity data the error 
bars are smaller than the symbol sizes. The horizontal lines show the solar reference abundance, 
from Anders & Grevesse (1989). 
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ABSTRACT 

Sightlines through high- and intermediate- velocity clouds allow measurements of ionic gas phase abun- 
dances, A, at very low values of H I column density, N(H I). Present observations cover over 4 orders of 
magnitude in N(H I). Remarkably, for several ions we find that the A vs N(H I) relation is the same at 
high and low column density and that the abundances have a relatively low dispersion (factors of 2-3) 
at any particular N(H I). Halo gas tends to have slightly higher values of A than disk gas at the same 
N(H I), suggesting that part of the dispersion may be attributed to the environment. We note that the 
dispersion is largest for Na I; using Na I as a predictor of N(H I) can lead to large errors. 

Important implications of the low dispersions regarding the physical nature of the ISM are: (a) because 
of clumping, over sufficiently long pathlengths N(H I) is a reasonable measure of the local density of most 
of the H atoms along the sight line; (6) the destruction of grains does not mainly take place in catastrophic 
events such as strong shocks, but is a continuous function of the mean density; (c) the cycling of the 
ions becoming attached to grains and being detached must be rapid, and the two rates must be roughly 
equal under a wide variety of conditions; (d) in gas that has a low average density the attachment should 
occur within denser concentrations. 

Subject headings: ISM: clouds, ISM, Abundances Galaxy: halo, Galaxy: disk 



1. INTRODUCTION 

The study of elemental abundances in the interstellar 
medium has a long history, both in the optical and with 
space-based UV spectrographs, such as Copernicus and 
the Hubble Space Telescope (HST). These studies show 
that most elements have apparent abundances relative to 
hydrogen substantially below those in the Sun. This is 
interpreted as due to depletion of ions onto dust grains. 

Savage & Sembach (1996a) summarized high-quality 
abundance studies based on HST data. They concluded 
that the depletion of an element depends on the nature of 
the interstellar medium along the sightline. In cold, dense 
disk gas, S has near-solar abundance, while Si, Mg, Mn, 
Cr, Fe and Ni are progressively more depleted. In warm- 
disk gas, the gas-phase abundances appear higher than in 
the cool gas, while in halo gas they are even larger. 

Jenkins, Savage & Spitzer (1986) used Copernicus to 
measure accurate column densities for Mg II, P II, CI I, 
CI II, Mn II and Fe II, for stars with distances, d, be- 
tween 0.1 and 2kpc. They measured N(H I) from the 
damping wings of Lya and found a relation between the 
ionic abundance relative to H I, A, and the average den- 
sity, (rin) (= N(H i)/d). An explanation was provided by 
Spitzer (1985), who suggested that the observed depletion 
in a line of sight depends on the fraction of warm and 
cold (dense) gas encountered, where each component has 
a specific gas phase abundance of each ion. 

Crinklaw, Federman, & Joseph (1994) observed Ca II 
and Ti II toward 12 stars with distances between 150 and 
650 pc. The slope of the relation between log[(riH)] and 
log[A(ion)] is similar for Ca II and Ti II, but it differs 
from those for Fe II, Mn II, Mg II and P II. They proposed 
that Ca II and Ti II mostly occur in the warm intercloud 



medium, and that the apparent depletion is very sensitive 
to the inclusion of high-density regions (with high deple- 
tion) in the sightline. 

These studies were limited to sightlines with N(H I) in 
the range 10 20 to 7xl0 21 cm -2 . To extend this range, we 
used the high- and intermediate-velocity gas (HVCs and 
IVCs, see Kuntz & Danly 1996, Wakker & van Woerden 
1997, Wakker 2000). In 21-cm emission these stand out 
in velocity from gas in the Milky Way disk, so that weak 
components (down to 10 18 cm -2 ) can be measured. At 
low velocities, parcels of gas with such low column densi- 
ties can not be separated from other high column density 
concentrations in the line of sight. 

Wakker (2000) analyzed all the published absorption- 
line data pertaining to the high- and intermediate- velocity 
gas. Improved H I column densities (and thus improved 
ion abundances) were determined for about 250 sightlines. 
Ca II was measured for 97 components, Ti II for 6, Fe II for 
13, Mn II for 5, Mg II for 10, and Na I for 34 (low-velocity 
Na I was taken from Ferlet et al. 1985). Most IVCs have in- 
trinsically near-solar abundance (Wakker 2000). However, 
HVC complex C has A~0A solar (Wakker et al. 1999), 
and the Magellanic Stream has A~0.25 solar (Gibson et 
al. 2000). Including these HVC complexes could lower 
some of the apparent abundances (in 8 cases for Ca II, 
4 for Fe II and 10 for Mg II), but no clear effect is visi- 
ble. [An anomalously low value for Mg II, in the VHVC 
probed by Mrk205, was excluded as this cloud may have 
low intrinsic abundance, and N(H I) is uncertain.] 

No distances are known for many of the HVCs/IVCs, 
and they fill only a small fraction of the pathlength. Thus, 
we looked at the relation between A and N(H I), rather 
than (nn)- Unexpectedly, we found tight anticorrelations 
of N(H I) with the abundances of Mg II, Ca II, Ti II, Mn II 
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and Fe II. A minimal correlation is expected because both 
quantities are integrated along lines of sight with varying 
physical conditions. Variations include (a) the ionization 
of H, (b) widely differing histories of grain destruction by 
interstellar shocks, and (c) considerable differences of the 
interstellar radiation field impinging on the gas. Ca II and 
Na I are not even the dominant stage of ionization in neu- 
tral regions, so the radiation field is of critical importance 
for their ionization fraction. In Sect. 2 we show the corre- 
lations between column density and abundances, while in 
Sect. 3 we discuss some implications. 

2. CORRELATIONS BETWEEN ABUNDANCES AND N(H I) 
2.1. The data 

Figure 1 shows scatter plots of log[N(H I)] vs log[A(ion)], 
the gas-phase abundance relative to H I, for Mg II, Ca II, 
Ti II, Mn II, Fe II and Na I. The least-squares fit and dis- 
persion (labeled "rms") are given in each panel. Symbols 
indicate the origin of the datapoints, as detailed in the 
figure caption. The Local ISM points (open squares) were 
not used in the fits (see Sect. 2.5 below). The high-quality 
data obtained with the "Goddard High Resolution Spec- 
trograph" (GHRS) are presented separately, in order to 
show the difference between halo and disk sightlines more 
clearly (Sect. 2.5). 

The currently preferred oscillator strengths of Mn II- 
1197, 1201 and Mg 11-1239, 1240 are 0.20 and 0.67 dex 
higher than those used by Jenkins et al. (1986) (see Sav- 
age & Sembach 1996a). We therefore corrected their Mn II 
and Mg II column densities downward by these amounts. 

As the highest abundances are associated with high- 
velocity gas, we might be seeing a manifestation of the 
Routly & Spitzer (1952) effect (this is the increase in 
the ratio N(Ca II)/N(Na I) with LSR velocity). This ef- 
fect is interpreted as showing that Ca is less depleted at 
higher peculiar velocities. However, it is based on nearby 
(<100pc), low- velocity (<20kms _1 ) gas (e.g. Vallerga et 
al. 1993). For the HVCs/IVCs the LSR velocity is not 
a good measure of the peculiar velocity relative to their 
surroundings. Further, the IVCs for which both Na I and 
Ca II have been measured do not show the effect (Wakker 
2000). 

We now discuss some implications of the A vs N(H I) 
relations. 

2.2. Result 1: small dispersions 

For Mg II, Mn II and Fe II the standard deviations of 
log[A(ion)] from the least-squares fit at a given N(H I) are 
only 0.27 in the log. That is less than a factor 2 either 
way, even though the correlation extends over 4 orders of 
magnitude in N(H I) and over 2 orders of magnitude in 
A(Fe II). The standard deviations are slightly larger for 
Ca II and Ti II (0.4 in the log, or a factor 2.6 either way). 
The largest scatter is seen for Na I, where it is 0.52 in the 
log or a factor 3.5 either way. 

The larger scatter for Ca II and especially Na I may be 
related to the fact that, unlike the other ions, these are not 
the dominant ionization stage in the diffuse ISM, so that 
ionization effects may play a larger role. However, Ti II is 
the dominant ionization stage, and it shows a scatter that 
is comparable to that of Ca II. 



Ferlet et al. (1985) proposed a slope near zero for the 
relation between Na I abundance and N(H I) for the low- 
velocity gas. This has been used extensively to estimate 
N(H I) from N(Na I). The lower column density IVC data 
tend to show slightly higher abundances. However, the 
spread in A(Na I) at any given value of N(H I) is rather 
large. For a given N(Na I), the derived N(H I) will be 
a factor >3.5 off either way in 33% of the cases. The 
most deviant points in the Na I diagram differ from the 
mean relation by a factor >20. Thus, although on average 
N(Na I)/N(H I) is independent of N(H I), N(Na I) is a 
rather poor predictor of N(H I). 

2.3. Result 2: differing slopes 

The slopes of the correlations differ substantially: 
-0.78±0.04 for Ca II, -0.69±0.08 for Ti II, -0.59±0.04 
for Fe II, -0.39±0.04 for Mn II, -0.24±0.04 for Mg II, and 
-0.16±0.06 for Na I. We find slopes of -0.63±0.12 and 
— 0.66±0.23 for low- and high-velocity Ca II; for Fe II we 
find -0.54±0.05 and -0.70±0.07, respectively. For other 
elements there are insufficient datapoints to make separate 
fits. The errors are larger because of the reduced range in 
N(H I), but the resulting slopes are the same to within the 
formal error. This shows that on average the high column 
density, low-velocity disk gas behaves in the same manner 
as the low column density high-velocity halo gas. 

2.4. Result 3: no obvious ionization effects 

The ionization potentials (LP.) of Fe II, Mn II, and 
Mg II are 16.2, 15.6, and 15.0 eV, respectively (compared 
to 13.6 eV for H I), so these ions can co-exist with H + . 
The Fe III/Fe II ratio should not depend on N(H I), but 
only on the surrounding radiation field (which only de- 
pends on location). Thus, if the interstellar radiation field 
has similar values in the sightlines that were observed, and 
if H+ were only present in shells around H I cores, then 
the fraction of H + should increase with decreasing N(H I), 
resulting in an overestimate of the abundance, when cal- 
culated as N(ion)/N(H I). This would produce an upturn 
at low N(H I) in the A vs N(H I) relation. However, such 
an effect is not obvious. Call (LP. 11.9eV) should ionize 
more easily than H, leading to a downturn at low N(H I). 
This is also not seen. 

2.5. Result 4-' environmental dependence 

For Mg II, Mn II and Fe II the GHRS disk points 
(open circles in Fig. 1) scatter around those of Jenkins 
et al. (1986). GHRS halo sightlines (filled circles) tend 
to lie above the average correlation. The small scatter in 
A(Fe II) remarked upon by Savage & Sembach (1996a) is 
consistent with the fact that their halo points span only 
a factor 10 in N(H I). However, some more recent halo 
points (filled triangles) have higher Fe Il/H I ratios, as do 
some HVCs. 

The low N(H I) Local ISM (LISM; open squares) shows 
abundances far below halo points at the same N(H I), al- 
though they are near the high end for disk gas and the ob- 
served depletion is consistent with the average Local ISM 
density of 0.1 cm -3 . We surmise that both the Sun and the 
LISM probes lie within a cloud, so that the LISM paths 
do not pass through all of it, unlike the case for clouds 
in the sightlines to more distant stars. Thus, the Local 
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ISM points suggest that density is the relevant physical 
parameter determining depletion. Over sufficiently long 
pathlengths column density can substitute because it is 
dominated by the densest regions. 

The largest depletions are seen toward ( Oph and £ 
Per, which were used by Sembach & Savage (1996) as the 
archetypical cold disk sightlines. These have high N(H I) 
and A(Mg II) and A(Fe II) lie below the average relation. 
Including N(H 2 ) would move these points 0.2 and 0.4 dex 
to the right and down, respectively, shifting them ~0.1 
and 0.2 dex further below the average. 

We thus see the environmental dependence found by 
Sembach & Savage (1996): disk sightlines tend to have 
low gas-phase abundances, while halo sightlines show high 
gas-phase abundances. We can understand this in large 
part as a column density dependence, as disk sightlines 
tend to have higher column density than halo sightlines. 

3. DISCUSSION 

The general trend that the gas phase abundances of 
refractory ions decreases with N(H I) is, of course, well 
known; the ions are removed from the gas phase by sticking 
to grains. However, the low dispersions of the depletions 
at a particular value of N(H I), continuing to low values to 
N(H I), are remarkable. The depletion patterns provide a 
diagnostic of the density fluctuation within the ISM, since 
the depletions require encounters between the grains and 
ions. The low dispersions of A at a given N(H I) suggest 
additional properties of the ISM: 

1. Large values of N(H I) arise from one or a few regions 
of relatively large local density rather than from a superpo- 
sition of many regions of lower density. In this way, there 
is unexpected physical meaning to the integrated column 
density: a good correlation with the mean local density 
of the typical H atom along the sightline. If large values 
of N(H I) were built up from many regions of low density 
along the line of sight, the averaged A would correspond 
to the lower densities. The association of large N(H I) 
with a large local density is strengthened by the fact that 
the anticorrelation of gas phase abundances with (nn) is 
better than with N(H I) (Savage & Bohlin 1979). 

2. The low dispersions in log[^4] suggest that grains are 
not destroyed in strong shocks, as some models have sug- 
gested. If such destruction occurred, a few points with 



large gas phase abundances would be expected at large 
N(H I). Rather, the detachment of a refractory element 
from grains seems to be a continuous function of the mean 
density, and is common enough to produce a low disper- 
sion about the mean. Grain destruction may occur by 
means of many small shocks with a broad distribution of 
velocities. 

3. Both the rates of ions sticking to grains and becom- 
ing detached from them (probably in atomic form) must 
be rapid enough to achieve rather good statistical equilib- 
rium. These balancing rates lead to a particular value of 
A for a given density; otherwise, the past history of the 
H I along the line of sight would strongly influence A. 

4. In view of (3) above, we think that the best possibility 
is that attachment takes place within regions of far greater 
than average density distributed along each line of sight. 
Encounters between grains and ions are very slow if the 
mean density of the diffuse ISM is assumed (see below). 
Possibly magnetohydrodynamic waves create non-thermal 
motions of the ions relative to the grains, since the charge 
per mass of the grains is vastly less that for the ions. As 
a result, ions encounter more dust grains and the net at- 
tachment rate increases. 

Under the assumption that relative motions were ther- 
mal, Weingartner & Draine (1999) calculated average gas 
phase abundances of a highly depleted species (Ti II) 
for large N(H I). In their model, dense ((riH)~10 3 cm -3 ) 
molecular clouds contain half the mass of the ISM, 30% 
is in cool, neutral clouds ((riH)=30 cm~ 3 ) and 20% is in 
warm neutral clouds ((riH)=0.4cm~ 3 ). The attachment 
of the Ti onto grains occurs in the molecular and cool 
clouds. This model shows that the observed A(Ti II) for 
large N(H I) can only be achieved by circulation into very 
dense regions. To achieve rapid attachment within the dif- 
fuse ISM, such circulation seems necessary, operating at 
both low and high column densities. An additional condi- 
tion for rapid attachment of ions is that the gas be cold as 
well as dense, in order to provide efficient focusing by the 
Coulomb force onto negatively charged grains. 

We appreciate comments and clarifications by an anony- 
mous referee, Alex Lazarian, Chris Howk, and Bruce 
Draine. B.P.W. is partially supported by NASA grant 
NAG5-9179. 
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Figure 1. Correlation plots of the abundances of Na I, Mg II, Ca II, Ti II, Mn II and Fe II vs N(H I). The solid lines show the 
least-squares fits, whose coefficients are given by the label in the lower left corner, "rms" is the rms of the residual, also indicated by 
the two parallel dotted lines, "p" is the correlation coefficient. Filled circles are for the HVC and IVC data from Wakker (2000). If 
measurement errors were given for N(ion), error bars are shown; usually no errors are available. Crosses show data for low- velocity 
gas from Ferlet et al. (1985; Na I), Jenkins et al. (1986; Mg II, Mn II and Fe II) and Crinklaw et al. (1994; Ca II, Ti II). Errors 
are usually smaller than the size of the dot. The right panels for Mg II, Mn II and Fe II separately show measurements with 
relatively low systematic errors obtained with the GHRS by Savage et al. (1993), Sofia et al. (1994), Cardelli et al. (1995), Spitzer 
& Fitzpatrick (1995), Savage & Sembach (1996b) and Sembach & Savage (1996). For these data, halo sightlines are shown by 
filled circles, disk sightlines by open circles. Filled/open triangles show halo/disk sightlines from Fitzpatrick & Spitzer (1997; after 
combining their fit components into the 5 H I components visible in the spectrum) and Howk et al. (1999). Open squares show 
Local ISM sightlines from Linsky et al. (1995), Hebrard et al. (1999) and Howk (priv. comm.); these were not used in the fits, see 
Sect. 2.4. For all low-velocity data the error bars are smaller than the symbol sizes. The horizontal lines show the solar reference 
abundance, from Anders & Grevesse (1989). 



